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ABSTRACT 

We report phase-referencing VLBA observations of H 2 0 masers near the star- 
forming region W3(OH) to measure their parallax and absolute proper motions. The 
measured annual parallax is 0.489 ± 0.017 milli-arcseconds (2.04 ± 0.07 kpc), where 
the error is dominated by a systematic atmospheric contribution. This distance is con¬ 
sistent with photometric distances from previous observations and with the distance 
determined from CH 3 OH maser astrometry presented in a related paper. We also find 
that the source driving the H 2 0 outflow, the “TW-object”, moves with a 3-dimensional 
velocity of > 7 km s^ 1 relative to the ultracompact H II region W3(OH). 
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1. Introduction 

The annual parallax is the most direct mea¬ 
surement of distances in astronomy. The Hippar- 
cos satellite successfully measured the distances 
to numerous stars in the Solar neighborhood, 
typically achieving 10 % accuracies for distances 
of ~ 100 pc, which contributed significantly to 
many fields of modern astronomy (e.g. Perry¬ 
man et al. 1995). However, annual parallax mea¬ 
surements for stars with kpc distances require 
sub-milliarcsecond accuracy, which has not been 
achieved optically. 

Very Long Baseline Interferometry (VLBI) 
provides the highest resolution in astronomy. In 
phase-referencing VLBI, the position of a target 
source is measured relative to a nearby positional 
reference source (see e.g. Beasley & Conway 
1995; Ros 2003). The feasibility of annual par¬ 
allax measurements with the Very Long Baseline 
Array (VLBA) has been demonstrated at low fre¬ 
quencies by Brisken et al. (2002) who measured 
annual parallaxes of pulsars in the Galaxy and 
by van Langevelde et al. (2000) and Vlemmings 
et al. ( 2002 ) who measured distances of Galac¬ 
tic OH masers associated with late type stars. 
Chatterjee et al. (2004) measured pulsar paral¬ 
laxes at 5 GHz and showed that the accuracy of 
astrometric measurements improves with higher 
frequencies. Their results indicate that one can 
measure distances of up to a few kpc with better 
than 10% uncertainty with VLBA astrometry of 
maser sources. Indeed, Kurayama et al. (2005) 
used the VLBA to measure the annual parallax 
of the Mira-Type star UX Cygni with high accu¬ 
racy. 

Hence, VLBA measurements allow sources spread 
over a large part of the Milky Way to have accu¬ 
rate parallaxes. This enables us to probe Galac¬ 
tic structure and dynamics since maser sources 
are spread over the whole Galaxy and, especially 
water vapor (H 2 O) maser sources are even found 
in its outer reaches (e.g. Wouterloot et al. 1993). 

The 22.2 GHz transition of H 2 O is the most 


widespread and luminous known maser line. In 
our Galaxy it has been detected toward numer¬ 
ous evolved red giant stars and high- and low- 
mass star-forming regions (see, e.g., Valdettaro 
et al. 2001 ). 

W3(OH) is a region containing several high- 
and intermediate-mass young stars and proto¬ 
stars of different evolutionary stages (e.g. Wilner 
et al. 1999 ; Wyrowski et al. 1997, 1999). 
In addition to strong OH and CH 3 OH masers, 
which are seen projected on the archetypical ul¬ 
tracompact H II (UCH II) region, very strong 
H 2 O maser emission is found toward the Turner- 
Welch (TW) Object (Turner & Welch 1984; Reid 
et al. 1995; Wilner et al. 1999), a protostar 
projected ~ 10 4 AU east of the UCH II region. 
The W3(OH) H 2 O masers were amongst the first 
studied with VLBI (Moran et al. 1973). VLBI 
Maps of the H 2 O maser emission have been re¬ 
ported by Alcolea et al (1992). 

We observed W3(OH) to measure its annual 
parallax and to study the internal dynamics of 
the known bipolar H 2 O outflow from the TW 
object. Moreover, our observations constitute 
a trial parallax and proper motion observation 
to explore the potential of utilizing H 2 O masers 
as probes of Galactic structure. Here we report 
VLBA observations of the W3(OH) H 2 O masers 
which yielded an extremely accurate parallax. 

2. Observations and data reduction 

We observed the W3(OH) H 2 O masers seven 
times with the NRAO VLBA 10 between January 
2001 and May 2002 (see Table 1). Each obser¬ 
vation was carried out over a 4-hour period in¬ 
cluding calibrator observations. The separations 
between the epochs were between two and four 
months. We observed two 16-MHz bands with 

10 The VLBA is operated by the National Radio Astron¬ 
omy Observatory (NRAO). The National Radio Astron¬ 
omy Observatory is a facility of the National Science Foun¬ 
dation operated under cooperative agreement by Associ¬ 
ated Universities, Inc. 
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one band centered on the maser velocity. The 
data were correlated with 1024 spectral channels 
in each band with an integration time of 2 sec¬ 
onds. The resulting velocity channel spacing was 
0.224 km s^ 1 , and we covered a velocity range of 
229 km s^ 1 . 


Table 1: Summary of the VLBA observations. 



Epoch 

Stations 

2001/01/28 

01:11:00 - 05:13:00 (UT) 

7 

2001/05/12 

18:20:00 - 22:20:00 (UT) 

9 

2001/07/12 

14:20:00 - 18:20:00 (UT) 

8 

2001/08/25 

11:27:00 - 15:27:00 (UT) 

9 

2001/10/23 

07:35:00 - 11:35:00 (UT) 

9 

2002/01/12 

02:16:00 - 06:16:00 (UT) 

9 

2002/05/06 

18:44:00 - 22:44:00 (UT) 

7 


We used ICRF 0244+624 as a phase-reference 
source. Its angular separation from W3(OH) is 
2.2°. Since this source is extragalactic with a 
redshift of 0.0438 (Margon & Kwitter 1978), its 
proper motion should be negligible. The source 
was detected at all epochs with peak flux den¬ 
sities > 0.9 Jy. The source was compact and 
unresolved as in previous VLBI observations at 
lower frequencies (Fey & Chariot. 2000), mak¬ 
ing it an excellent phase-reference source. Typ¬ 
ical source elevations varied from 46 to 62 de¬ 
grees. We switched every 20 seconds between the 
sources W3(OH) and ICRF 0244+624, yielding 
typical on-source times of 7 seconds. The strong 
source NRAO150 was observed for 5 minutes ev¬ 
ery 44 minutes for delay and bandpass calibra¬ 
tion. 

The data were calibrated and imaged with 
standard techniques using the NRAO Astronom¬ 
ical Image Processing System (AIPS) software 
package. Amplitude calibration used system 
temperature measurements and standard gain 
curves. A fringe fit was performed on ICRF 
0244+624 and the solutions were applied to W3(OH). 
The Kitt Peak and Los Alamos antennas did not 
observe in the first epoch because of heavy snow, 
the Pie Town antenna was flagged in the third 
and seventh epoch since most of the data was lost 


because of system failures. Also, the Saint Croix 
antenna was flagged in all epochs as it produced 
little useful data, probably owing to its moist, 
low-latitude, low-elevation site and the absence 
of short baselines to this telescope. 

Given an accurate geometric model in the 
VLBA correlator, the largest errors in phase- 
referencing observations are introduced by a zenith 
delay error in the atmospheric model of the corre¬ 
lator (see Reid et al. 1999). These errors degrade 
the image quality and the astrometric accuracy. 
The main contributions to the fringe-phase of the 
target source, after phase-referencing, are from a 
position offset and the atmospheric/ionospheric 
delay error (if the source structure is negligible). 
Because of the different behavior of the two con¬ 
tributions, it is possible to separate both effects 
and to estimate the position offset as well as a 
zenith delay error. We fit a model phase owing 
to a position offset and a zenith delay error at 
each station to the calibrated phase data of an 
isolated and strong maser feature. The phase er¬ 
rors caused by the zenith delay errors can then be 
corrected by the AIPS task CLCOR. This correc¬ 
tion improves the quality of the phase-referenced 
images and the astrometric accuracy (Reid et al. 
1999; Reid & Brunthaler 2004; Brunthaler et al. 
2005). The data from the Mauna Kea antenna 
at the 7th epoch was flagged, since we could not 
accurately estimate the zenith delay error of this 
station. After applying these corrections to the 
UV data of W3(OH) we created images of the 
maser features for each epoch. 

We determined the position of each maser spot 
by fitting a two-dimensional Gaussian brightness 
distribution to the maps, using the AIPS task 
JMFIT. The formal error of the positions were 
calculated from the signal-to-noise ratio of the 
fitted peak flux densities and sizes of the compo¬ 
nents and were typically 10 /was in right ascension 
and 20 +as in declination. 
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3. Results 


We found a total of 42 distinct H 2 O maser 
features that were detectable over 3 epochs, usu¬ 
ally in 2 to 10 adjacent velocity channels with 
peak flux densities from a few hundred mJy to 
a few hundred Jy (Tab. 2). The masers were 
distributed over an area of 2.5” x 0.5” (Fig. 1), 
which is consistent with previous VLBI observa¬ 
tions (Alcolea et al. 1992). 

The absolute proper motion of a maser feature 
is the sum of the outflow motion in W3(OH), 
annual parallax, Galactic rotation, Solar motion 
and peculiar motion of W3(OH). We assumed 
all motions except the annual parallax to be lin¬ 
ear. All motions, except the internal motions, are 
equal for all maser features. It can be challeng¬ 
ing to trace exactly the same maser feature, since 
H 2 O masers are highly time variable and their 
absolute proper motions relative to the extra- 
galactic reference source are non-linear because 
of the effect of the parallax. For the feature 
identification we used not the absolute proper 
motions but the motions relative to a reference 
feature (feature 1 in Table 2) and then looked 
for rectilinear motions of a reasonable magnitude 
(i.e., < 100 km s _1 ) in each spectral channel. 
The relative motions in Table 2 were obtained 
after a phase self-calibration on the reference fea¬ 
ture to reduce systematic errors. 

3.1. The annual parallax of W3(OH) 

Once we identified a maser feature over five or 
more epochs, we modeled its path in terms of its 
proper motion (^ a , /is) and the annual parallax 
(II) by using 

Aacosd = Hf a (a, 5,t) + Li a t + ao 
A5 = nf s (a,5,t) + nst + 5 0 

where t is time, ao and 5q are the positions of 
a maser feature at t = 0, and the functions f a 
and fs are the parallax displacements in right 
ascension and declination, respectively, given, for 
example, by Smart (1965). 



East Offset fmasl 

Fig. 1.— H 2 O maser distribution. Filled cir¬ 
cles show features detected in 7 epochs, open cir¬ 
cles features detected in 6 or 5 epochs detectable, 
and open triangles show features detected in 4 or 
3 epoch. The numbers denote the components 
listed in Table 1 for which astrometrical fits were 
obtained. 



Fig. 2.— Spectra of component 7 at epoch 2 
(exes), epoch 3 (stars), epoch 4 (open squares), 
epoch 5 (filled squares), epoch 6 (open circles), 
and epoch 7 (filled circles). 
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Table 2 

Detected H 2 0 masers 


No. 

Offset 3 

E-W N-S 

(mas) 

epochs 

Relative proper motion* 3 
E-W N-S 

(mas yr -1 ) 

V LSR C 

Max Min 
(km s -1 ) 

Peak flux d 
Max Min 

(Jy/beam) 

1 

-144.49 

-144.29 

7 

0 

±0 

0 

±0 

-50.9 

-51.3 

304.6 

57.7 

2 

2.33 

6.27 

6 

0.62 

0.13 

2.30 

0.18 

-47.9 

-48.8 

95.8 

12.4 

3 

1.97 

6.52 

5 

0.54 

0.28 

2.63 

0.19 

-47.9 

-48.2 

165.5 

33.3 

4 

-4.92 

-5.68 

7 

1.06 

0.06 

1.26 

0.09 

-48.2 

-48.8 

335.5 

65.4 

5 

-9.45 

3.00 

4 

0.97 

0.14 

1.41 

0.19 

-48.2 

-48.6 

253.1 

21.6 

6 

-9.63 

33.42 

7 

0.37 

0.06 

2.73 

0.05 

-48.4 

-49.0 

367.0 

40.5 

7 

-8.22 

33.69 

6 

0.87 

0.11 

2.57 

0.11 

-48.6 

-49.6 

127.0 

12.6 

8 

-6.38 

33.18 

5 

-0.15 

0.11 

3.06 

0.07 

-49.8 

-49.8 

12.9 

5.5 

9 

-4.83 

32.42 

4 

0.66 

0.58 

2.88 

0.18 

-49.6 

-49.8 

20.6 

4.9 

10 

-19.26 

39.29 

4 

-0.17 

0.17 

2.81 

0.13 

-48.8 

-49.4 

72.3 

9.9 

11 

-13.67 

49.22 

5 

0.83 

0.08 

1.51 

0.20 

-49.2 

-49.6 

255.2 

21.9 

12 

-21.64 

63.22 

6 

1.09 

0.29 

1.02 

0.18 

-49.0 

-49.2 

126.1 

51.9 

13 

-20.30 

62.46 

5 

0.42 

0.15 

2.16 

0.09 

-48.8 

-49.0 

170.1 

49.1 

14 

-21.95 

64.34 

4 

0.40 

0.31 

2.20 

0.62 

-49.0 

-49.2 

77.8 

18.1 

15 

-26.84 

75.60 

5 

0.47 

0.07 

1.84 

0.20 

-48.4 

-48.4 

55.2 

9.7 

16 

-38.68 

98.35 

6 

0.01 

0.09 

2.01 

0.04 

-49.6 

-49.6 

102.4 

14.5 

17 

-643.08 

90.22 

6 

0.08 

0.23 

1.49 

0.08 

-51.1 

-51.5 

41.8 

00 

bo 

18 

-648.91 

91.76 

3 

2.95 

0.10 

1.82 

0.13 

-53.0 

-53.2 

9.7 

5.7 

19 

-650.25 

92.27 

5 

1.86 

0.34 

1.94 

0.39 

-52.6 

-51.7 

10.4 

3.7 

20 

-650.64 

92.11 

5 

0.73 

0.10 

1.43 

0.10 

-52.6 

-52.8 

35.7 

3.7 

21 

-663.46 

95.49 

6 

0.47 

0.08 

1.40 

0.06 

-54.3 

-54.7 

22.5 

3.2 

22 

-674.24 

96.97 

3 

3.96 

0.27 

1.11 

0.14 

-55.1 

-55.1 

7.7 

3.1 

23 

-673.21 

96.43 

3 

4.59 

0.30 

1.04 

0.41 

-54.9 

-55.1 

2.1 

1.4 

24 

-680.99 

95.14 

4 

0.07 

0.09 

1.51 

0.10 

-55.7 

-56.0 

6.6 

1.0 

25 

-758.57 

-48.82 

3 

6.86 

0.07 

-3.12 

0.15 

-59.5 

-60.2 

3.5 

0.6 

26 

-826.70 

-38.99 

3 

3.40 

0.04 

-0.80 

0.01 

-57.2 

-57.2 

7.2 

2.6 

27 

-827.14 

-38.23 

5 

3.75 

0.14 

-0.53 

0.06 

-57.4 

-57.6 

2.6 

0.9 

28 

-831.22 

-35.39 

3 

1.26 

1.04 

-0.70 

0.01 

-56.8 

-57.4 

3.2 

1.6 

29 

-863.41 

-12.12 

5 

2.58 

0.23 

-0.03 

0.12 

-58.5 

-60.2 

30.0 

1.0 

30 

-942.98 

1.95 

3 

-4.60 

0.34 

-0.18 

0.71 

-64.6 

-65.0 

0.6 

0.2 

31 

-963.36 

59.78 

3 

-4.01 

0.99 

4.00 

0.15 

-51.9 

-52.4 

32.9 

18.6 

32 

-964.18 

60.34 

3 

-3.28 

0.24 

3.07 

0.64 

-51.7 

-51.9 

57.2 

12.9 

33 

-975.72 

54.48 

5 

-5.02 

0.14 

3.01 

0.13 

-59.1 

-60.8 

2.0 

0.8 

34 

-976.17 

54.40 

6 

-4.77 

0.37 

3.31 

0.10 

-59.5 

-61.2 

2.6 

0.9 

35 

-982.77 

80.99 

3 

-1.44 

0.14 

1.19 

0.10 

-53.6 

-54.1 

283.8 

121.9 

36 

-983.46 

81.71 

3 

-0.06 

0.06 

-3.85 

0.11 

-53.4 

-53.6 

658.0 

230.5 

37 

-983.47 

82.56 

4 

-1.14 

0.08 

0.81 

0.33 

-53.0 

-53.4 

70.5 

9.3 

38 

-985.41 

-24.01 

3 

-3.86 

0.04 

0.82 

0.07 

-53.8 

-54.1 

97.7 

52.0 

39 

-988.70 

75.89 

4 

-2.91 

0.09 

2.81 

0.08 

-61.4 

-61.8 

3.5 

0.6 

40 

-1034.87 

-34.44 

3 

-3.20 

0.05 

0.26 

0.13 

-57.8 

-58.1 

6.5 

2.1 

41 

-2045.95 

-133.41 

7 

-1.93 

0.06 

1.36 

0.09 

-45.2 

-45.2 

28.7 

7.8 

42 

-2231.40 

84.12 

6 

-3.48 

0.16 

0.67 

0.17 

-62.5 

-63.1 

121.1 

4.3 


a Origin of position offsets is a J2 ooo = 2 h 27 m 04.8362±0.0004 s , <5j 20 oo = 61° 52 24.607i0.002". 

b Motion of strongest channel relative to the reference feature 1. Errors are normalized for \ 2 ■ 

5 

c Time variation of the LSR velocity in the channel with the highest peak flux. The velocity drift 
is constant for some H 2 0 masers, while it is fluctuating within listed velocities for others. This is 
most likely caused by blending with the intensities of blended features changing relative to each 
other. Full line width for each feature is from 0.6 to 3.4 km s“ . 



First, we fitted a proper motion and an an¬ 
nual parallax to all velocity channels individu¬ 
ally. The reduced x 2 values of the fits were very 
high (10 - 20). This was caused by unrealistically 
small formal errors of the position estimates, es¬ 
pecially for strong sources. Hence, we introduced 
an error floor by adding quadratically a value of 
0.05 mas to the formal position error. Possible 
sources of the error floor are variation of the cen¬ 
troid position of extragalactic source, residual er¬ 
rors in the estimation of zenith delay corrections 
and blending of maser features. This resulted a 
reduced x 2 near unity. Thus, the positional ac¬ 
curacy of a single channel in a single epoch is ~ 
50 fias. 

It is very difficult to quantify the individual er¬ 
ror contributions. The variations in the centroid 
position of the extragalactic reference source can 
be caused by unresolved structure changes and 
could be a few tens of microarcseconds per year. 
This motion should not influence the parallax 
measurement, if the motion is constant over the 
time our experiment. However, in the case of 
an ejection of a new jet component one could 
get non-linear motions that affect the parallax 
measurements. This can not be excluded since 
we did not use a second extragalactic reference 
source in our observations. However, ejections of 
new jet components are usually accompanied by 
sharp rises in the flux density. We find only small 
variations of the flux density of ICRF 0244+624 
(~ 10%) and consider this scenario as unlikely. 

Residual errors in the estimates of zenith delay 
corrections are much more likely. Since our ob¬ 
servations were relatively short (~ 4 hours) and 
do not cover a large range of different elevations, 
it is difficult to separate the atmospheric and 
position offset contributions in the fringe phase. 
These residual errors can lead to position errors 
of several tens of microarcseconds for an indi¬ 
vidual observation. Since the atmospheric con¬ 
ditions between epochs are not correlated the re¬ 
sulting position errors are random. 

A major problem in the parallax fitting is that 


most maser features show strong variability be¬ 
tween the observations. Indeed, most features 
were not detected in all seven epochs. The flux 
densities of the maser features often change by 
Tip to a factor of 10, typically causing significant 
changes in the (blended) line shape and shifts 
in the apparent center velocity (Fig. 2). Since 
there is strong variability in blended spectra, one 
expects some effect on the astrometric accuracy 
of our measurements. Indeed, the fits yielded 
a spread in parallaxes with typical values be¬ 
tween 0.45 and 0.55 mas, and a few outliers at 0.4 
and 0.6 mas. This scatter in the parallax values 
is much larger than the formal accuracy. This 
scatter can not be explained by the two previ¬ 
ously mentioned sources of error (centroid posi¬ 
tion variation of the reference source and atmo¬ 
spheric contributions), since they affect all maser 
components similar. 

For some components, the parallax model could 
not match the measured positions, as seen in re¬ 
duced x 2 values larger than ~ 3. (In these cases, 
we probably are dealing with physical changes 
of the source, such as a brightening of one part 
of its only marginally resolved structure. Also, 
a new component with a similar radial veloc¬ 
ity might appear at nearly the same position, 
while the existing one fades.) Other components 
had good fits, as evidenced by reduced x 2 val¬ 
ues near unity, but still showed a large scatter in 
the parallaxes. This can be explained by a cor¬ 
relation between the parallax and proper motion 
parameters, coupled with position errors caused 
by structural changes in some masers. Signifi¬ 
cant correlation coefficients are the result of non- 
optimal time sampling of position measurements, 
caused by telescope scheduling and maser vari¬ 
ability. The effect of the non-optimal time sam¬ 
pling is not identical for all maser components, 
because not all components were detected in all 
epochs. 

The fits of most features show large correlation 
coefficients between the parallax and the proper 
motions. A histogram of the multiple correlation 
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Table 3: Best fit annual parallax for H 2 O masers in W3(0H). The errors of the annual parallaxes are 
from 20 to 50 //as for each spot, while the global fit shows 9 //as statistically. 


Comp. 

v LS r [km s" 1 ] 

7 r [mas] 

1 

-50.3 

0.516 ± 0.037 

1 

-50.5 

0.479 ± 0.025 

1 

-50.7 

0.489 ± 0.021 

1 

-50.9 

0.494 ± 0.020 

1 

-51.1 

0.497 ± 0.020 

1 

-51.3 

0.490 ± 0.020 

1 

-51.5 

0.482 ± 0.018 

Average 

0.492 ± 0.009 

4 

-48.2 

0.465 ± 0.047 

4 

-48.4 

0.486 ± 0.038 

Average 

0.476 ± 0.029 

Global fit 

0.489 ± 0.009 



Multiple correlation coefficient 


Fig. 3.— Histogram of the multiple correlation 
coefficients of the parallax parameter for differ¬ 
ent maser features. The dashed horizontal line 
marks the upper limit for the coefficient used to 
determine the final parallax. 


coefficient of the parallax clearly shows a bimodal 
distribution (Fig. 3) - with a few values below 0.3 
and a large number with values larger than 0.4. 
Hence, we used only components that were de¬ 
tected in all seven epochs and showed a multiple 
correlation coefficient of the parallax parameter 
of < 0.3. This left seven channels of component 
1 , which is the strongest and spectrally broadest 
of all features, and two channels of component 
4. These components also had very symmetric 
spatial brightness distributions. The spectra of 
the two features in all seven epochs are shown 
in Fig. 4. The parallaxes of the individual chan¬ 
nel fits are in good agreement within their joint 
errors (Table 3). 

We also obtained a global fit to all nine chan¬ 
nels with one parallax, but allowing a different 
proper motion for each channel. The measured 
positions and our model for two channels are 
plotted in Fig. 5. This global fit yields a parallax 
of 0.489 ± 0.009 mas. The 0.009 mas uncertainty 
is statistical only and does not include systematic 
effects. Xu et al. (2005) see indications that par¬ 
allax measurements might show some systematic 
sensitivity to the angular offset of the calibra- 
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Fig. 4.— Spectra of component 1 (upper panel) 
and 4 (lower panel). See caption of Fig. 3 for 
details. 


tors. For this effect, they included a systematic 
parallax error of about 0.007 mas per degree of 
calibrator separation in their total uncertainty. 
This is probably caused by residual errors in the 
zenith delay correction. To be conservative, since 
ICRF 0244+624 has a separation from W3(OH) 
of 2.2 degrees, we add a systematic component 
of 0.015 mas to the statistical uncertainty of 
0.009 mas. Therefore, we find the annual par¬ 
allax of W3(OH)-TW to be 0.489 ± 0.009 ± 0.015 
mas where the first error indicates the statistical 
error from the measurements while the second 
error describes a systematic atmospheric contri¬ 
bution. 

To investigate this systematic error in more 
detail, we performed simulations where we cal¬ 
culated the position offsets given a parallax (0.5 
mas), random proper motions (between -4 and 
+4 mas yr” 1 ), position of W3(OH) on the sky, 
and dates of the observations. Then we added 
a random Gaussian error with an rms of 50 //as 
(the position accuracy of a single channel in a 
single epoch in our data) to each position offset, 
and fitted the simulated data set. We fitted 1000 
simulated data sets and found that the result¬ 
ing parallaxes followed a Gaussian distribution 
with a standard deviation of ~ 0.015 mas. These 
errors also affect the proper motion fits. The dif¬ 
ference between the true proper motion and the 
fitted proper motions in our simulations has a 
standard deviation of 0.04 mas yr -1 . 

The annual parallax for W3(OH)-TW cor¬ 
responds to a distance of 2.04+ 0.04 ± 0.06 
kpc. This is far more accurate than any pre¬ 
vious distance, and comparable to the result of 
Xu et al (2005). Photometric distance estimates 
of 2.2 kpc to an OB association near W3(OH) 
(Humphreys 1978) compare favorably with our 
value. However, the kinematic distance for a 
source at Galactic longitude 133.95 degrees with 
an LSR velocity near —50 km s _1 is 5 kpc. The 
reason for this large discrepancy is a peculiar mo¬ 
tion of the W3(OH) region, which is discussed in 
detail in Xu et al. (2005). 




3.2. Relative maser motions in W3(OH)- 
TW 


The relative proper motions of the different 
maser features can be described by an outflow 
model. This has been successfully applied to 
earlier VLBI observations of the H 2 O masers in 
W3(OH) by Alcolea et al. (1993, hereafter A93). 
To estimate the physical parameters of the out¬ 
flow model of W3(OH), we used the proper mo¬ 
tions, radial velocities and positions of maser fea¬ 
tures (Table 2).We chose component 1 from Ta¬ 
ble 2 as a reference feature and used its data 
to re-calibrate the all of the maser data. Com¬ 
ponent 1 was chosen because it is one of the 
strongest and least affected by blending of the 
masers (see Sect. 3.1). We used only the maser 
features that were detected in at least 3 epochs. 
Fig. 6 (lower panel) shows the relative internal 
proper motions. We used the same method as 
A93, but we adopted our accurate distance of 
2.04 kpc. We estimated the velocity and posi¬ 
tion of the center of expansion with respect to 
the reference maser feature and the expansion 
velocity V exp at 1” from the center of expansion. 
Details of the model fitting are described in A93 
and Irnai et al. (2000). The best fit to the data 
was obtained minimizing the expression 



where (u x , u y , u z ) are the motions of the maser 
feature in right ascension, declination and the 
radial velocity, W x and W y are the tangential 
motions relative to the reference feature and W z 
is the radial velocity of the center of expansion. 
(v x ,v y ,v z ) are components of a maser spot’s ve¬ 
locity, v, which is given by the equation of v 
= Vexp \r\ a r/ \r\ where r = (x - X 0 , y - Y 0 , z). 
(a x , (Jy, (J z ) are components of the measurement 
uncertainty and a possible turbulent velocity of 
4 km s -1 (added in quadrature). Finally Xq and 


Yo are the positions of the center of expansion rel¬ 
ative to the phase tracking center. We assigned 
a value of —51.0 km s _1 for the LSR velocity of 
the center of expansion (W 2 ), based on thermal 
molecular emissions from the location of the TW 
object (Wyrowski et al. 1997). 

The results of the best model fit with a re¬ 
duced x 2 °f 2.7 are shown in Table 4. The for¬ 
mal errors were increased by y/2.7 to account for 
the x 2 . It is not possible to compare the po¬ 
sition of center of expansion between A93 and 
our results, since the A93 results were not phase- 
referenced. However, the parameter a agrees 
well with the value of A93. We find a value of 
13.0 ± 5 km s' 1 for the expansion velocity which 
is slightly lower than the 20 ± 2 km s -1 value 
obtained by A93. This discrepancy is probably 
not significant, since A93 used a larger distance 
of 2.2 kpc that leads to larger velocities. When 
scaling the A93 results to our distance, one gets 
an expansion velocity of 18.5 ± 2 km s -1 . Then 
the difference is 5.5 ± 5.4 km s -1 . 

In our case, the absolute position of the ref¬ 
erence maser feature can be determined with re¬ 
spect to the extragalactic source, and transfered 
to the other masers features and the center of 
expansion. Fig. 6 displays the water masers and 
their tangential motions on the 8.4 GHz Very 
Large Array (VLA) continuum map of Wilner 
et al. (1999) and the 220 GHz Plateau de 
Bure interferometer (PdBI) continuum map of 
Wyrowski et al. (1999). Absolute positions of 
the TW object and the center of expansion of the 


Table 4: Best fit model for H 2 O maser velocity 


field 


Parameter 

X 0 

-0.89 

± 

0.05 

arcsec 

Y 0 

0.00 

± 

0.03 

arcsec 

w* 

-10.0 

± 

7 

km s _1 

Wy 

16.0 

± 

7 

km s _1 

V exp 

13.0 

± 

5 

km s _1 

a 

-0.38 

± 

0.2 
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H 2 O maser outflow are listed in Table 5. The 
absolute position of center of expansion is con¬ 
sistent with the TW object, in which the H 2 O 
maser outflow originates. 

3.3. Relative motion with respect to the 
Ultracompact H II region in W3(OH) 

The tangential motion of the reference feature 
relative to the extragalactic background source 
ICRF 0244+624 obtained from the parallax/proper 
motion fit is —1.2 ± 0.5 km s _1 eastward and 
— 10.2 ± 0.5 km s' 1 northward, for a distance 
of 2.04 kpc. Adding the tangential motion of 
the center of expansion relative to these values 
gives the tangential motion of the center of ex¬ 
pansion: — 11.2 ± 7 ± 0.4 km s^ 1 eastward and 
5.8 +7+0.4 km s^ 1 northward. The second error 
of 0.4 km s' 1 comes from the systenratics dis¬ 
cussed in Sect. 3.1. While this systematic error 
dominates the error in the parallax measurement, 
the fitting uncertainty of the outflow model dom¬ 
inated the total error in the proper motion of the 
center of expansion. 

On the other hand, similar astrometric phase- 
referencing VLBA observations have been car¬ 
ried out for 12.2 GHz methanol masers associ¬ 
ated with the UCH II region in W3(OH) (Xu 
et al. 2005). These authors have estimated the 
proper motion of these methanol masers with re¬ 
spect to extragalactic continuum sources, and ob¬ 
tained a tangential motion of —11.1 ± 0.2 km s^ 1 
and —1.3 ± 0.1 km s -1 , eastward and northward 
respectively. They have not estimated the inter¬ 
nal motions of the methanol masers. However 
these are approximately 2 km s _1 (Moscadelli et 
al. 2002), much smaller than the internal H 2 O 
maser velocities. Thus we assign a total uncer¬ 
tainty of 2 km s^ 1 to the absolute tangential mo¬ 
tion of the methanol masers. 

Combining the two VLBA results, the cen¬ 
ter of expansion of the H 2 O masers (presum¬ 
ably the TW object) moves at —0.1 ±7.3 km s^ 1 
eastward, +7.1 ± 7.3 km s^ 1 northward and —7 
km s -1 in radial direction, with respect to the 


UCH II region of W3(OH). We used a systemic 
radial velocity of -51 km s -1 for the H 2 O masers 
(TW object) and -44 km s _1 for the UCH II 
region in W3(OH). Converting this relative mo¬ 
tion to a Galactic Cartesian coordinate system, 
the TW object relative to the UCH II region 
moves toward the Galactic center with a veloc¬ 
ity of 7.3 ± 5.9 km s -1 , in the direction opposed 
to Galactic rotation with 2.1 ± 4.3 km s -1 , and 
toward the North Galactic Pole with 6.5 ± 7.3 
km s -1 . 

4. Discussion 

4.1. 3D dynamics and structure of W3(OH) 

In the case of W3(OH), a total (3-dimensional) 
motion of the TW object with respect to the 
UCH II region is > 7 km s -1 (7 km s _1 in radial 
velocity plus 7+10 km s -1 tangential motion). 
What causes this large relative motion? Assum¬ 
ing that the TW object and the UCH II region 
are gravitationally bound, a total mass (Aft) of 
the W3(OH) region can be estimated by 


> 1.1 x 10 3 ( / ) ( - ———pM 0 , 

V10 4 AU / V10 km s 1 / 

where r is a separation of W3(OH)-TW and 
UCH II region and v is the relative motion. Since 
the separation along the line of sight is still un¬ 
known, a separation of 10 4 AU is a minimum 
value. This total mass is much higher than any 
reasonable estimate for the combined mass of 
the stars in W3(OH)-TW the UCH II region, 
for which we estimate 17 and 23 Mg, respec¬ 
tively, since the spectral type of the TW object 
is approximately B0 (Wyrowski et al 1999) and 
that of the UCH II is approximately 08.5 region 
(Harten 1976). So, these objects do not appear 
to be gravitationally bound. 

The only other massive stars with similarly 
measured 3-dimensional motions are the Becklin- 
Neugebauer object and radio source-I (related 


10 



Table 5: Absolute positions of the continuum peak and the HoO maser center of expansion. 


Position 

R.A. (J2000) 

Dec.(J2000) 

8.4 GHz 

02 27 04.713 

61 52 24.65 

220 GHz 

02 27 04.71 

61 52 24.6 

Center of expansion 

02 27 04.7103 ±0.0071 

61 52 24.607T0.030 


For the peak positions in the 8.4 and 220 GHz images we assume uncertainties of 0.02 and 0.1 arcseconds, typical for high 
quality VLA and PdBI images. See text for the determination of the H 2 O center of expansion uncertainty. 
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Fig. 5.— Positions versus time for two maser 
components (Vlsr of -50.9 km s -1 for compo¬ 
nent 1 and -48.4 km s -1 for component 4) rela¬ 
tive to the compact extragalactic source ICRF 
0244+624. The top (bottom) plots show the 
northward (eastward) position components. Best 
fit parallax and proper motions are indicated 
with dashed lines. 
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to IRc 2) in the Orion Kleinmann-Low region. 
The relative motion of these objects is very large, 
> 45 km s -1 (Rodriguez et al. 2005), and Tan 
(2004) invokes ejection of the Becklin-Neugebauer 
object from the core of the Orion Nebula Trapez¬ 
ium Cluster to explain this. However, for the 
W3(OH) sources, given their present-day rela¬ 
tive motions, a close encounter of TW and the 
UCH II region cannot have happened. Perhaps, 
one or the other (or both) had an encounter with 
a third object in the recent past and are now un¬ 
bound. 

4.2. Limitations of H 2 O maser astrome¬ 
try 

The accuracy of the present H 2 O maser paral¬ 
lax measurements achieved with phase-referencing 
VLBI is ~ 10 /xas. However, the high time vari¬ 
ability of H 2 O masers limits such measurements 
in significant ways. While, we detected over 40 
maser features at any epoch, and were able to 
trace 20 features over at least 5 epochs, only 2 of 
these 20 maser features had measurements that 
yield a reliable parallax measurement. In typical 
sources, few features might persist over a period 
of > 1 yr, which is optimum for annual parallax 
measurements. 

If one wishes to use the parallax and proper 
motion results to study Galactic structure and 
kinematics, one needs to model the internal mo¬ 
tions of the masers. To obtain an accurate model 
fit of the internal motions one needs to mea¬ 
sure the motions of many maser features. Poor 
estimation of the internal motion (typical mo¬ 
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Fig. 6.— Upper panel: the 8.4 GHz VLA continuum map of W3(OH) shown in contours (Wilner et al. 
1999) with the absolute tangential motions of TW object and the W3(OH) UCH II region indicated with 
cones. The colors indicate radial velocities as coded in the panel at the lower right. The gray cone shows 
their relative tangential motion. Lower panel: the internal tangential motions of the W3(OH)-TW H 2 O 
masers. The origin of the lower figure is (aj 2 ooo = 2 /l 27 m 04.8362 s , <5j2ooo = 61° 52' 24.607"). The thin 
contours and grey scale represent the 220 GHz Plateau de Bure interferometer continuum map (Wyrowski 
et al. 1999). The peak positions for the TW object in this and the 8.4 GHz map are consistent with the 
expansion center of the outflow model. 
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tions are 20 to 200 km s^ 1 ) leads to inaccurate 
estimates of the 3-dimensional motion. Of the 
many hundreds of known H 2 O maser sources in 
the Milky Way, most will not have as many de¬ 
tectable components as W3(OH). Thus it may be 
difficult to study Galactic dynamics using only 
H 2 O masers. 

Due to strong variability and large internal ve¬ 
locities, H 2 O maser sources are not the best can¬ 
didates to study the Galactic structure and dy¬ 
namics. Other maser sources that show less vari¬ 
ability and slower internal motions (e.g. methanol 
maser Xu et al. 2005) are preferable. On the 
other hand, H 2 O masers are much more com¬ 
mon and also found in Galactic regions where no 
methanol masers are found. Some H 2 O masers 
are found in the outer Galaxy near the edge of 
the optical stellar disk (e.g. Wouterloot et al. 
1988), while methanol masers (e.g. Pestalozzi et 
al. 2005) have not been found there. Therefore, 
these H 2 O masers can be the best sources to mea¬ 
sure distances and motions in the outer Galaxy. 

5. Conclusions 

We have measured the annual parallax of the 
H 2 O maser source in the W3(OH) region with 
phase-referenced VLB A observations. The dis¬ 
tance of 2.04 ± 0.07 kpc that we obtain is con¬ 
sistent with previous photometric distance esti¬ 
mates (but with much higher accuracy) and with 
the CH 3 OH maser parallax corresponding to 1.95 
± 0.04 kpc determined by Xu et al. (2005) in the 
related paper. 

We also measured the proper motions of the 
W3(OH)-TW H 2 O masers and find that the TW 
object is moving with a speed of > 7 km s _1 with 
respect to the nearby UCH II region (with its OH 
and CH 3 OH masers). Such a large speed differ¬ 
ence between two massive objects in the same 
star forming region is puzzling. 

Although H 2 O masers are not perfect target 
sources to investigate Galactic structure and dy¬ 
namics, they can still provide important infor¬ 


mation about regions in the Galaxy that are not 
accessible otherwise (e.g. the outer Galaxy). 
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